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We study the statics and dynamics of a dipolar Bose-Einstein condensate (BEC) droplet bound 
by inter-species contact interaction in a trapped non-dipolar BEC. Our findings are demonstrated 
in terms of stability plots of a dipolar ^^^Dy droplet bound in a trapped non- dipolar ^'^Rb BEC with 
a variable number of ^^^Dy atoms and the inter-species scattering length. A trapped non-dipolar 
BEC of a fixed number of atoms can only bind a dipolar droplet containing atoms less than a 
critical number for the inter-species scattering length between two critical values. The shape and 
size (statics) as well as the small breathing oscillation (dynamics) of the dipolar BEC droplet are 
studied using the numerical and variational solutions of a mean-field model. We also suggest an 
experimental procedure for achieving such a ^^^Dy droplet by relaxing the trap on the ^^^Dy BEC 
in a trapped binary ^^Rb-^^^Dy mixture. 
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I. INTRODUCTION 

The experimental observation of a dipolar Bose- 
Einstein condensate (BEC) of ^^Cr ^^^Dy 
and ^^^Er [8] atoms with large magnetic dipole moments 
has opened new directions of research in cold atoms 
in the quest of novel and interesting features related 
to the anisotropic long-range dipolar interaction. Polar 
molecules, with much larger (electric) dipole moment, 
are also being considered [9|] for BEC experiments. The 
atomic interaction in a dilute BEC of alkali-metal and 
other types of atoms (with negligible dipole moment) is 
represented by an S-wave contact (delta-function) poten- 
tial. However, the non-local anisotropic long-range dipo- 
lar interaction acts in all partial waves and is attractive 
in certain directions and repulsive in others. 

An untrapped three-dimensional (3D) BEC with at- 
tractive interaction does not exist in nature due to col- 
lapse instability [l^. However, the collapse of an un- 
trapped BEC can be avoided due to interspecies contact 
interaction in a binary mixture with a trapped BEC. The 
shape of such a stable droplet bound in a trapped BEC 
is controlled by the inter-species interactions, whereas 
that of a trapped BEC is determined by the underlying 
trap. Here we consider such a dipolar droplet bound in 
a trapped nondipolar BEC. The effect of the underly- 
ing atomic interaction, specially that of the anisotropic 
dipolar interaction, will easily manifest in such a dipolar 
droplet. These droplets will be termed quasi- free as they 
can easily move around inside the larger trapped BEC re- 
sponsible for their binding. By taking the trapped BEC 
to be nondipolar, one can easily study the effect of intra- 
species dipolar interaction on the quasi-free droplet in 
the absence of any inter-species dipolar interaction. 

Dipolar BECs are immediately distinguishable from 



those with purely contact interactions by their strong 
shape and stability dependence on trapping geome- 
try. Here, we are proposing a new way to trap a dipolar 
BEC in the form of a quasi- free dipolar droplet, using an 
attractive inter-species mean-field potential, which intro- 
duces a unique trapping geometry that results in further 
interesting shape and stability characteristics of dipolar 
BECs. We study the statics and dynamics of a quasi- 
free dipolar droplet using numerical solution and varia- 
tional approximation of a mean-field model BEl. For 
this purpose we consider a binary mixture of non-dipolar 
Rb87 and dipolar Dyl64 where the trapped non-dipolar 
^^Rb BEC could be in a cigar- or disk-shape. Among 
the available dipolar BECs |3, 0, Q, "^^"^Dy atoms have 
the largest (magnetic) dipolar interaction strength ^ . 
The existence of stable 3D ^^^Dy droplets is illustrated 
by stability plots involving the number of atoms of the 
two species and the inter-species scattering length. For a 
fixed number A^(Rb) of ^^Rb atoms the dipolar droplet 
could be bound below a critical number iV(Dy) of ^^^Dy 
atoms between two Hmiting values of inter-species attrac- 
tion. If the inter-species attraction is too small, the ^^^Dy 
atoms in the droplet cannot be bound and, for a very 
large inter-species attraction, the droplet is destroyed by 
collapse instability. Usually, the dipolar droplet has the 
same shape as the trap acting on the non-dipolar BEC. 
However, for a large number of ^^^Dy atoms, as the inter- 
species attraction approaches the collapse instability, the 
dipolar droplet always is of cigar shape even if the exter- 
nal trap is disk shaped and eventually the dipolar droplet 
collapses on the axial z axis from the cylindrical configu- 
ration. For a small number of atoms, the dipolar droplet 
collapses to the center maintaining its shape rather than 
on the axial z axis. 
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The variational approximation to the sizes and chemi- 
cal potentials of the stationary droplets is compared with 
the numerical solution of the mean-field model. The 
numerical study of breathing oscillation of the stable 
dipolar droplet is found to be in reasonable agreement 
'users/adhikari f with a time-dependent variational model calculation. We 
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also demonstrate a viable experimental way of creating 
a ^^^Dy droplet bound in a trapped ^^Rb BEC, e.g., by 
slowly removing the trap on a trapped binary ^^^Dy-^^Rb 
mixture, while the trapped ^^^Dy BEC evolves into a 
quasi- free droplet. 

In Sec. II the mean-field model for a dipolar BEC 
droplet bound in a trapped non-dipolar BEC is devel- 
oped. A time-dependent, analytic, Euler-Lagrange Gaus- 
sian variational approximation of the model is also pre- 
sented. The results of numerical calculation are shown in 
Sec. III. Finally, in Sec. IV we present a brief summary 
of our findings. 



II. MEAN-FIELD MODEL FOR A QUASI-FREE 
DIPOLAR DROPLET 



We consider a binary BEC, where one of the species is 
dipolar and the other non-dipolar, interacting via inter- 
and intra-species interactions with the mass, number, 
magnetic /colorred dipole moment, and scattering length 
for the two species i = 1,2, denoted by rrii^Ni^ jii^ai^^ 
respectively. The first species (^^Rb) is taken to be non- 
dipolar (/ii = 0) and trapped while the second species 
(^^^Dy) is dipolar (/i2 7^ 0) and polarized along axial 
z direction. The angular frequencies for the axially- 
symmetric trap on ^^Rb along x, y and z directions are 
taken as ujx = oOy = uoi and Uz = Xiooi. The inter- 
(V12) and intra-species (Vi) interactions for two atoms at 
positions r and are taken as 

Vi2(R) = R , ^i(R = -din), (1 
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where R = r — r', /io is the permeability of free space, 6 
is the angle made by the vector R with the polarization 
z direction, and mu = 17111712/ {mi + 777,2) is the reduced 
mass of the two species of atoms. With these interac- 
tions, the coupled Gross-Pitaevskii (CP) equations for 
the binary dipolar BEC can be written as |13| 
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To compare the dipolar and contact interactions, the 
intra-species dipolar interaction is expressed in terms of 
the length scale add, defined by add = /ioM2^2/(127rfi^). 
We express the strength of the dipolar interaction in Eq. 
(jlj by this length scale and transform Eqs. ([3]) and Q 
into the following dimensionless form [l3| 
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where mi2 = mi/m2, gi = 47raiA^i, g2 = 47ra2A^2^i2, 
gi2 = 27rmiai2N2/mR, g2i = 27rmiai2Ni/mR, gdd = 
3A^2<^dd^i2- In Eqs. (|6]) and ([7j), length is expressed 
in units o f oscillator length of the first species Iq = 
y/h/miUJi^ energy in units of oscillator energy huJi^ den- 
sity in units of and time in units of to = cjf ^. 

Convenient analytic variational approximation to Eqs. 
([6]) and ([7j) can be obtained with the following ansatz for 
the wave functions [l#jl6l 
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where r = {p, = {x,?/}, Wpi and Wzi are the widths 
and ai and Pi are additional variational parameters. The 
effective Lagrangian for the binary system is 
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where A = w^^ + w'^p2^ ^ = '^^i + ^ = 

4ai2mi/(y^mi^). In these equations the overhead dot 
denotes time derivative. The Euler-Lagrange variational 
equations for the widths for the effective Lagrangian (p!Q|) , 
obtained in usual fashion fl6j, are: 
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The solution of the time-dependent Eqs. ([TT]) — ([T3|) 
gives the dynamics of the variational approximation. 

If \ii is the chemical potential with which the stationary 
wave function ^^(r,t) propagates in time, e.g. (j)i{r^t) ~ 
ex.p{—ifiit)(j)i{T), then the variational estimate for /i^ is: 
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The widths of the stationary state can be obtained from 
the solution of Eqs. (pTj) — (p!4|) setting the time deriva- 
tives of the widths equal to zero. This procedure is equiv- 
alent to a minimization of the energy ([19]), provided the 
stationary state is stable and corresponds to a energy 
minimum. 



III. NUMERICAL RESULTS 

We solve Eqs. ([6]) and ([7]) by split-step Crank-Nicolson 
discretization scheme using a space step of 0.1 and the 
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FIG. 1: (Color online) |ai2| — A^(Dy) stability plot showing 
the domain of stable ^^^Dy droplet in a binary ^^Rb-^^^Dy 
mixture for (a) N{Rh) = 2000, (b) N(Rh) = 10000, and (c) 
A^(Rb) = 50000, and for A = 0.25 and 4. The intra-species 
scattering lengths are taken as ai = llOao- 



time step 0.001 [HI, The contribution of the dipolar 
interaction is calculated in momentum space by Fourier 
transformation [15]. For both species of atoms ^^Rb and 
^^^Dy we take the intra-species scattering length as = 
llOao, and the strength of dipolar interaction as add = 
131ao 0, The yet unknown inter-species scattering 
length ai2 is taken as a variable. The variation of ai2 can 
be achieved experimentally by the Feshbach resonance 
technique [18]. We consider the trap frequency ui = 
27r X 115 Hz, so that the length scale Iq = ^Jh/mlUOl = 1 
jim and time scale = uj^^ = 1.38 ms. 

We find that a quasi-free ^^^Dy droplet is achievable for 
a moderately attractive inter-species attraction (negative 
ai2) and for appropriate values of the number of atoms 
of the two species A^i and N2. We illustrate the domain 
of existence of a stable ^^^Dy droplet in terms of stability 
plots in Figs. [11(a), (b) and (c) for 7V(Rb) = 2000, 10000 
and 50000, and for A = 4 (disk-shaped trap) and 0.25 
(cigar-shaped trap). We consider cigar- and disk-shaped 
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FIG. 2: (Color online) Variational (v) and numerical (n) re- 
sults for chemical potential /x and rms sizes (x), (z) for the 
trapped ^'^Rb BEC of 10000 atoms and the bound ^^^Dy 
droplet versus the number of ^^'^Dy atoms for (a) A = 4 
and (b) A = 0.25. The inter-species scattering length is 
ai2 = — 30ao. 



i^^Dy droplets in this study, where the effect of dipolar 
interaction is expected to be more prominent, and not 
spherically symmetric ones where the effect of dipolar 
interaction is expected to be a minimum [H, [5[. In all 
plots of Figs. [H for a fixed N{Rh), the ^^^^Dy droplet 
can be bound for a maximum of N(Dy). This maximum 
of N(Dy) increases with A^(Rb), all other parameters re- 
maining fixed. For a fixed A/'(Rb), the disk-shaped trap 
can accommodate a larger maximum number of ^^^Dy 
atoms than the cigar-shaped trap. For A^(Dy) smaller 
than this maximum number, the ^^^Dy droplet can be 
bound for |ai2| between two limiting values. For |ai2| 
above the upper limit, there is too much inter-species at- 
traction on the the ^^^Dy droplet leading to its collapse. 
For |ai2| below the lower limit, there is not enough at- 
traction and the ^^^Dy droplet cannot be bound. The 
lower limit of |ai2| is small and tends to zero as A^(Rb) 
tends to infinity. For fixed values of 7V(Rb) and 7V(Dy) 
and for small |ai2| near the lower limit, the disk-shaped 
configuration is favored and it can bind the ^^^Dy droplet 
for smaller values of |ai2|. However, for a fixed 7V(Rb) 
and for a large |ai2| near the upper limit, the disk-shaped 
configuration is favored only for a large iV(Dy) allowing a 
i^^Dy droplet for larger |ai2|, whereas the cigar-shaped 
configuration is favored for a small 7V(Dy) allowing a 
i^^Dy droplet for larger |ai2|. 

In Figs. [2] we illustrate the numerical and varia- 
tional results for chemical potentials and root-mean- 
square (rms) sizes {x) and (z) of the trapped ^^Rb BEC 
of 10000 atoms and of the bound ^^^Dy droplet versus 



TABLE I: Root-mean-square sizes of the binary ^^Rb-^^^Dy 
system of 50000 ^^Rb atoms from variational approximation 
()lip — (|14p (var), and an approximate energy minimization 
(approx). 
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Ar(Dy) for (a) A = 4 and (b) 0.25. Considering that the 
i^^Dy droplet may not have a Gaussian shape as assumed 
in the variational approximation, the agreement between 
the numerical and variational results is quite satisfactory. 

The existence of the quasi-free ^^^Dy droplet can be 
studied qualitatively by a minimization of the energy 
(p^ . However, the widths of the ^^Rb BEC for a fixed 
A^(Rb) do not vary much as A/'(Dy) or a\2 is varied. 
Hence for a qualitative understanding of the existence 
of a stable ^^^Dy droplet for a fixed 7V(Rb), we can make 
further approximation in the energy (p!9|) and take the 
widths of the trapped ^^Rb BEC to be the same as the 
widths of the ^^Rb BEC in the absence of ^^^Dy atoms 
under otherwise identical conditions. The widths of the 
i^^Dy droplet so obtained for 7V(Rb) = 50000 are com- 
pared in Table I for several values of ai2, A^(Dy), and 
A with the widths obtained from exact energy minimiza- 
tion. We find from Table I that the sizes of the ^^Rb for 
a fixed 7V(Rb) remains reasonably constant with respect 
to the variation of A^(Dy) and a\2 and that the approx- 
imate energy minimization with fixed widths for ^^Rb 
BEC leads to a good approximation to the widths of the 
i^^Dy droplet. In the numerical solution of the binary 
CP equations (|6]) and ([7j) we also verified that the shape 
and size of the trapped ^^Rb BEC are fairly independent 
of the presence or absence of the ^^^Dy droplet. Hence 
in the study of the shape, size and dynamics of the dipo- 
lar ^^"^Dy droplet bound in the trapped ^'^Rb BEC, the 
trapped BEC will only have a passive role and we shall 
highlight only the shape, size and dynamics of the ^^^Dy 
droplet in the following. 

We consider the shape of a stable ^^^Dy droplet in 
the stability plots of Figs. [1] and study its change as 
the collapse instability is approached. In Fig. [3] (a) we 
show the isodensity contour plot (density |(/)i(r)p) of a 
disk-shaped (A = 4) ^'^Rb BEC of 10000 atoms and that 
(density |^2(r)P) of the bound ^^^Dy droplet in Figs. [3] 
(b) for 7V(Dy) = 1000 (ai2 = -30ao), (c) for 7V(Dy) 
= 1000 (ai2 = -SOao), and in (d) for A/'(Dy) = 4000 
(ai2 = -30ao). The isodensity of the ^'^Rb BEC IS prac- 
tically the same in all three cases. Of the isodensities of 



FIG. 3: (Color online) (a) The isodensity contour of a disk- 
shaped (A = 4) ^'^Rb BEC of 10000 atoms and that of the 
bound ^^'^Dy droplet of (b) 1000 atoms (ai2 = -30ao), (c) 
1000 atoms (ai2 = -80ao), (d) 4000 atoms (ai2 = -30ao). 
The isodensity contour of the ^^Rb BEC is practically the 
same in all three cases. The density on the contour in all 
cases is 0.001/q ^ and length is measured in units of lo (= 1 



FIG. 4: (Color online) (a) The isodensity contour of a cigar- 
shaped (A = 0.25) ^'^Rb BEC of 10000 atoms and that of the 
bound ^^^Dy droplet of (b) 1000 atoms (ai2 = -30ao), (c) 
1000 atoms (ai2 = — 65ao), (d) 3000 atoms (ai2 = — 30ao). 
The isodensity contour of the ^^Rb BEC is practically the 
same in all three cases. The density on the contour in all 
cases is 0.001/^^ and length is measured in units of lo (= 1 



the ^^^Dy droplet, the one in Fig. [3](b) is deep inside the 
stability region far away from collapse instability. The 
parameters of Figs. [3] (c) and (d) are close to the region 
of collapse instability. Of these two, Fig. [3] (c) corre- 
sponds to a medium number of ^^^Dy atoms and Fig. [3] 
(d) to a large number. Independent of the initial shape, 
a ^^^Dy droplet with a medium number of ^^^Dy atoms 
always collapses towards the center with shrinking size. 
However, a ^^^Dy droplet containing a large number of 
atoms, independent of the associated trap symmetry, al- 
ways first takes a cigar shape and then collapses on the 
axial z direction. A strong dipolar interaction prohibits 
a collapse to center of a large ^^^Dy droplet and favors 
a cigar shape. In Figs. [3] the trap is disk-shaped. The 
medium-sized ^^^Dy droplet of Fig. [3] (c) collapses to 
center maintaining the disk shape as the net dipolar in- 
teraction is smaller in this case. The large ^^^Dy droplet 
in Fig. [3] (d) has changed its shape from the supporting 
disk-shaped trap to a cigar shape. If the attraction is 
further increased by increasing |ai2|, the ^^^Dy droplet 
of Fig. [3] (d) would collapse from the cigar shape on the 
axial z axis. 

In Fig. m (a) we illustrate the isodensity contour of 
a cigar-shaped (A = 0.25) ^'^Rb BEC of 10000 atoms 
and the same of the ^^"^Dy droplet are shown in Figs. 
[1(b) for 7V(Dy) = 1000 (ai2 = -30ao), (c) for A^(Dy) 
= 1000 (ai2 = -65ao), and in (d) for N{Dy) = 3000 
(ai2 = -30ao). The density of the ^^Rb BEC is virtually 
the same in all three cases. Of the three ^^^Dy droplets, 
the one in Fig. [H (b) is deep inside the stability region far 
away from collapse instability. The parameters for Figs. 
[3] (c) and (d) are close to the region of collapse instability. 
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FIG. 5: (Color online) The rms sizes (a) (p) and (b) (r) 
from numerical simulation (num) and variational approxima- 
tion (var) during breathing oscillation initiated by the sudden 
change ai2 1.01 x ai2 versus time in units of to (= 1-38 
ms) of a ^^^Dy droplet of 100 atoms bound in a trapped ^^Rb 
BEC of 2000 atoms for ai2 = -50 and A = 4. 



In all cases the ^^^Dy droplet maintains the cigar shape 
of the accompanying trap acting on the ^^Rb BEC. The 
size of the ^^^Dy droplet in Fig. 2] (b) is larger than that 
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FIG. 6: (Color online) The rms sizes (a) (z) and (b) (x) 
from numerical simulation (num) and variational approxima- 
tion (var) during breathing oscillation initiated by the sudden 
change ai2 1.01 x ai2 versus time in units of to (= 1.38 ms) 
of a ^^^Dy droplet of 1000 atoms bound in a trapped ^^Rb 
BEC of 10000 atoms for ai2 = -50 and A = 0.25. 



of Fig. m (c) with the same number of atoms due to a 
strong inter-species attraction acting on the latter. The 
i^^Dy droplet of Fig. |4] (d) is larger than those in Figs, 
m (b) and (c) due to a larger number of ^^^Dy atoms in 
it. 

We investigate the dynamics of a bound ^^^Dy droplet 
in a trapped ^^Rb BEC. First, to test the present scheme 
we consider a small ^^^Dy droplet of 100 atoms bound 
in a disk-shaped (A = 4) ^'^Rb BEC of 2000 atoms with 
ai2 = — 50ao. This corresponds to the stable region of 
Fig. [T](a). In real-time evolution of the system, the scat- 
tering length ai2 is suddenly changed to 1.01 x ai2 thus 
starting the breathing oscillation. In Figs. [5]we show the 
resultant oscillation of the rms sizes (a) (p) and (b) (r) 
during time evolution as obtained from numerical solu- 
tion of Eqs. ^ and ([7j) and variational approximation 
(ED - (HI. 

After obtaining a satisfactory result of dynamics with 
a small system, we venture with a larger system of ex- 
perimental interest, e. g., a stable ^^"^Dy droplet of 1000 
atoms bound in a cigar-shaped (A = 0.25) ^''Rb BEC 
of 10000 atoms with ai2 = — 50ao corresponding to the 
stable region of Fig. [1] (b). In real-time evolution of 
the system, the scattering length ai2 is again suddenly 
changed to 1.01 x ai2 thus starting the breathing oscil- 
lation. In Figs. [6] we show the resultant oscillation of 
the rms sizes (a) (z) and (b) (x) during time evolution as 
obtained from a numerical solution of Eqs. (|6]) and ([7]) 



FIG. 7: (Color online) (a) The rms sizes (x) and (z) of the 
^^"^Dy and ^^Rb BECs versus time during the passage of the 
trapped ^^"^Dy BEC in the binary ^'^Rb-^^^Dy mixture with 
N{Rh) = 10000, A^(Dy) = 1000, ai2 = -50ao and A = 4 as 
the confining trap on the ^^^Dy BEC is relaxed for t > 10 ex- 
ponentially as Vtrap exp[— (t — 10)]Vtrap, SO that this trap is 
practically zero for t > 20. The linear ID densities (b) |(/?(2;)|^ 
and (c) \(p{x)\'^ at times t = 20, 40, 60, 80, 100 (in units of to) 
during the expansion together with the numerically calculated 
density of the stationary (sta) ^^'^Dy droplet. 



and variational approximation (pTj) — ([Hj). The agree- 
ment between the numerical simulation and variational 
approximation in both cases of breathing oscillation il- 
lustrated in Figs. [5] and [6] for disk- and cigar-shapes is 
quite satisfactory considering the facts that (i) during 
oscillation the bound ^^^Dy droplet might have a den- 
sity distribution different from a Gaussian distribution 
assumed in the variational approximation and (ii) the 
large non-linearity in the ^^Rb BEC will also make its 
density distribution non-Gaussian. The stable oscilla- 
tion under small perturbation as shown in Figs. [5] and [6] 
confirms the dynamical stability of the ^^^Dy droplet for 
both disk and cigar type environments. 

The present quasi- free dipolar droplet is not just of the- 
oretical interest, but can be realized experimentally by 
initially preparing a binary ^^Rb-^^^Dy mixture where 
both the components are harmonically trapped. The 
trap on ^^^Dy can then be ramped to zero exponen- 
tially during a few milli seconds. The ^^^Dy cloud will 
initially expand and finally form the quasi-free ^^^Dy 
droplet. To illustrate numerically the viability of this 
procedure, we consider an initial ^^Rb-^^^Dy binary mix- 
ture with following parameters: A/'(Rb) = 10000, A^(Dy) 
= 1000, ai2 = — 50ao. We take the same axial trap 
with A = 4 acting on both components. Then we per- 



7 



form real-time propagation of the binary GP equation 
with this initial state and ramp down the ^^^Dy trap by 
nrap(Dy) ^ exp[-(t - 10)]Vtrap(Dy) for time t > 10. 
The trap on ^^^Dy is practically zero for times t > 20. 
The trapped ^^^Dy BEG first expands for t > 10 and 
eventually it emerges as the quasi- free ^^^Dy droplet. 
The time evolution of the rms sizes of the ^^"^Dy droplet 
is shown in Fig. [3(a). The widths of the ^^^Dy droplet 
first increase and eventually execute small oscillation il- 
lustrating the stable ^^^Dy droplet. To investigate the 
shape of the ^^^Dy droplet, we plot in Figs. [71 (b) 
and (c) the one-dimensional (ID) densities |v^(^)P = 
j j dxdy\(l){x,y,z)\^ and \^{x)\^ = J J dydz\(J){x,y, z)\'^ 
at times t = 20, 40, 60, 80, and 100 along with the cor- 
responding numerical densities of the stationary droplet 
calculated from the solution of Eqs. ([6]) and (Q. The ID 
numerical densities of the stationary ^^^Dy droplet are 
in agreement with the densities obtained from dynamical 
simulation of the passage of the trapped ^^^Dy BEG to a 
quasi-free droplet. In Figs. [3 (b) and (c) one can easily 
identify the small oscillation of the evolving dynamical 
droplet around its stationary shape. 

IV. SUMMARY AND DISCUSSION 

Using variational approximation and numerical solu- 
tion of a set of coupled mean- field GP equations, we 
demonstrate the existence of a stable dipolar ^^^Dy 
droplet bound by inter-species attraction in a trapped 
non-dipolar BEG of ^^Rb atoms. The domain of stabil- 
ity of the ^^^Dy droplet is highlighted in stability plots 
of number of ^^^Dy atoms and inter-species scattering 
length ai2 for both cigar- and disk-shaped traps acting 
on the ^^Rb BEG. Results of variational approximation 
and numerical solution for the statics (sizes and chemical 
potentials) and dynamics (breathing oscillation) of the 
i^^Dy droplet are found to be in satisfactory agreement 
with each other. We also demonstrate numerically that 
such droplets can be obtained experimentally by consid- 
ering a trapped binary ^^Rb-^^^Dy BEG and then remov- 
ing the trap on ^^^Dy. The ^^^Dy BEG then expands and 
transforms into a bound dipolar droplet. 

Dipolar interaction among atoms is quite different from 
normal short-range atomic interaction and manifests in 
different ways in a trapped dipolar BEG. However, in a 
trapped dipolar BEG, the confining constraints could be 
too strong and could make the effect of dipolar inter- 



action difficult to observe. Special experimental set up 
and theoretical formulation might be necessary to study 
the effect of dipolar interaction. However, it will be much 
easier to see the effect of dipolar interaction in the present 
binary mixture of non-dipolar ^^Rb and dipolar ^^^Dy. 
Strong dipolar interaction in the axial polarization z di- 
rection should elongate the dipolar BEG along this direc- 
tion thus transforming it to a cigar shape. However, it 
will be difficult to observe this change in the presence of a 
harmonic trap along z. The effect of dipolar interaction is 
clearly seen in the present quasi- free ^^^Dy droplet of Fig. 
[3] (d) where due to a strong dipolar interaction the shape 
of the ^^^Dy droplet has changed from the disk shape — 
the shape of the ^'^Rb BEG of Fig. [3] (a) responsible for 
its binding — to the cigar shape. The thin disk-shaped 
^^Rb BEG stays only near the central z = plane of the 
^^^Dy droplet. Most of the ^^^Dy atoms lying outside 
the ^''Rb BEG are bound due to the intra-species dipolar 
interaction. In the dipolar ^^^Dy droplet of Fig. [3] (d) 
the dipolar interaction is playing a more important role 
in its binding and shape. 

The present quasi-free dipolar droplet may also find 
other interesting applications in the BEG phenomenol- 
ogy. Among the interesting features in a trapped dipolar 
BEG, one can mention its peculiar shape and stability 
properties ^l| , and D-wave collapse jl9| , anisotropic soli- 
ton, vortex soliton [g^ and vortex lattice [2i|, anisotropic 
shock and sound wave propagation [22| , anisotropic Lan- 
dau critical velocity [23|, stable checkerboard, stripe, 
and star configurations in a two-dimensional (2D) opti- 
cal lattice as stable Mott insulator [24] as well as super- 
fluid soliton [25I states. It would be of great interest to 
find out how these features and properties of a trapped 
dipolar BEG would manifest in a 3D quasi-free dipolar 
droplet. For example, a quasi- free dipolar droplet could 
be used in the experimental study of anisotropic sound 
and shock- wave propagation [22[ , collapse dynamics [l9| , 
anisotropic Landau critical velocity |23j , formation of 
vortex dipoles and vortex lattice [21| etc in a different 
setting of confinement which will facilitate the observa- 
tion of the effect of anisotropic dipolar interaction. 
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